Experimental Determination of In-Medium Cluster Binding Energies and Mott Points 

in Nuclear Matter 
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In medium binding energies and Mott points for d, t, "^He and a clusters in low density nuclear 
matter have been determined at specific combinations of temperature and density in low density 
nuclear matter produced in collisions of iTA MeV ''"Ar and ^''Zn projectiles with ^^^Sn and ^^''Sn 
target nuclei. The experimentally derived values of the in medium modified binding energies are in 
good agreement with recent theoretical predictions based upon the implementation of Pauli blocking 
effects in a quantum statistical approach. 
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I. INTRODUCTION 

Simple nuclear statistical equilibrium models assume 
that properties of the species in equilibrium are the same 
as those of the isolated species. While this assumption 
is tenable at very low density, it is untenable at higher 
densities where in-medium effects lead to dissolution of 
the clusters and a transition to cluster-free nuclear mat- 
ter. To deal with this intermediate density range, Typel 
et al. [H have developed a quantum statistical approach 
which includes cluster correlations in the medium and 
interpolates between the exact low-density limit and the 
very successful relativistic mean field (RMF) approaches 
appropriate near the saturation density. The general- 
ized RMF model developed attributes the decrease of 
the cluster fractions at high densities to a reduction of 
the cluster binding energies due to the Pauli blocking. 
This leads to the Mott effect of vanishing binding 
Well-defined clusters appear only for densities below ap- 
proximately 1/10 of the saturation density and get dis- 
solved at higher densities. The maximum cluster density 
is reached around the Mott density. Because of the pres- 
ence of strong correlations in the scattering state contin- 
uum that are effectively represented by one resonance, 
there is a non- vanishing cluster fraction above the Mott 
density [3] ■ We report here the first experimental deriva- 
tion of temperature and density dependent binding ener- 



gies of d, t, '^He and a clusters, directly from experimental 
particle yields. Experimental values for Mott points for 
d, t, ■^He and a clusters are in good agreement with the 
predictions made in reference 



EXPERIMENTAL TECHNIQUES 

We reported in Refs. ^ and \d\ that measurements of 
nucleon and light cluster emission from the participant 
matter which is produced in near Fermi energy heavy 
ion collisions could be employed to probe the EOS at 
low density and moderate temperatures where cluster- 
ing is important. The NIMROD Att multi-detector at 
Texas A&M University has now been used to extend our 
measurements to higher densities. Cluster production 
in collisions of 47A MeV '^•'Ar with "2,i24gj^ ^^^^j 642^ 
with "2,i24gj^ ^g^g studied. NIMROD consists of a 166 
segment charged particle array set inside a neutron ball 
6]. The charged particle array is arranged in 12 rings 
of Si-Csl telescopes or single Csl detectors concentric 
around the beam axis. The Csl detectors are 1-10 cm 
thick Tl doped crystals read by photomultiplier tubes. A 
pulse shape discrimination method is employed to iden- 
tify light particles in the Csl detectors. Neutron mul- 
tiplicity is measured with the An neutron detector sur- 
rounding the charged particle array. The combined neu- 
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tron and charged particle multiplicities were employed 
to select the most violent events for subsequent analysis. 
Further details on the detection system, energy calibra- 
tions and neutron ball efficiency may be found in refer- 
ence M. 



ANALYSIS 

The dynamics of the collision process allow us to probe 
the nature of the intermediate velocity "nucleon-nucleon" 
emission source Measurement of emission cross 

sections of nucleons and light clusters together with suit- 
able application of a coalescence ansatz Q provides the 
means to probe the properties and evolution of the inter- 
action region. The techniques used have been detailed 
in several previous publications 0, 0, 0- U | and are de- 
scribed briefly below. A notable difference from refer- 
ences and H is the method of density extraction. This 
is discussed more extensively in the following. We em- 
phasize that the event selection is on the more violent 
collisions. Cross section weighting favors mid-range im- 
pact parameters. 

An initial estimation of emission multiplicities at each 
stage of the reaction was made by fitting the observed 
light particle spectra assuming contributions from three 
sources, a projectile- like fragment (PLF) source, an inter- 
mediate velocity (IV) source, and a target-like fragment 
(TLF) source. A reasonable reproduction of the observed 
spectra is achieved. Except for the most forward detec- 
tor rings the data are dominated by particles associated 
with the IV and TLF sources. The IV source velocities 
are very close to 50% of the beam velocity as seen in 
many other studies (0-[l3l and references therein). The 
observed s pec tral slopes reflect the evolution dynamics of 
the source For further analysis, this IV source 

is most easily sampled at the intermediate angles where 
contributions from the other sources are minimized. For 
the analysis of the evolution of the source we have se- 
lected the data in ring 9 of the NIMROD detector. This 
ring covered an angular range in the laboratory of 38° to 
52°. The results of the three-source fit analyses, as well 
as inspection of invariant velocity plots constructed for 
each ejectile and each system, indicate that this selection 
of angular range minimizes contributions from secondary 
evaporative decay of projectile like or target like sources 

We treat the IV source as a nascent fireball created 
in the participant interaction zone. The expansion and 
cooling of this zone leads to a correlated evolution of 
density and temperature which we probe using particle 
and cluster observables, yield, energy and angle. As in 



the previous work 
tope yield ratios [14 



we have employed double iso- 
15j | to characterize the temperature 
at a particular emission time. Model studies comparing 
Albergo model temperatures and densities to the known 



input values have shown the double isotope ratio temper- 
atures to be relatively robust in this density range jl6| . 
However the densities extracted using the Albergo model 
are useful only at the very lowest densities Both of 
these results are confirmed in the more extensive calcu- 
lations of reference [H. In this study we have employed a 
different means of densityextraction, the thermal coales- 
cence model of Mekjian (sl. [lol|. 

To determine the coalescence parameter Pq, the radius 
in momentum space, from our data we have followed the 
Coulomb corrected coalescence model formalism of Awes 
et al. 17 1 and previously employed by us in reference Q ■ 
In the laboratory frame the derived relationship between 
the observed cluster and proton differential cross sections 
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where the double differential multiplicity for a cluster of 
mass number A containing Z protons and N neutrons 
and having a Coulomb-corrected energy Ea, is related 
to the proton double differential multiplicity at the same 
Coulomb corrected energy per nucleon, E — Ec, where 
Ec is the Coulomb barrier for proton emission. Rnp is 
the neutron to proton ratio. Since within the framework 
of the coalescence model the yield ratios of two isotopes 
which differ by one neutron are determined by their bind- 
ing energies and the n/p ratio in the coalescence volume, 
we have used the observed triton to "^He yield ratio to 
derive the n/p ratio used in this analysis. 

In the Mekjian model thermal and chemical equilib- 
rium determines coalescence yields of all species. Under 
these assumptions there is a direct relationship between 
the derived radius in momentum space and the volume 
of the emitting system. In terms of the Pq derived from 
Eq. (1) and assuming a spherical source 
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where h is Plancks constant and Z, N, and A are the 
same as in Eq. (1), Eq is the binding energy and s the 
spin of the emitted cluster and T is the temperature. 
Thus the volume can be derived from the observed Pq 
and temperature values assuming a spherical shape in 
terms of the Pq derived from Eq. (1). 

Because our goal was to derive information on the den- 
sity and temperature evolution of the emitting system, 
our analysis was not limited to determining an average 
Pq value. Instead, as in our previous studies 
results for d, t, '^He, and ^He, were derived as a function 
of Vsurf. the velocity of the emerging particle at the nu- 
clear surface, prior to Coulomb acceleration fv^ . From 
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the relevant Pq values we then determined volumes using 
Eq. (2). A comparison of these volumes indicated good 
agreement for t, ^He and ''He. The volumes derived from 
the deuteron data are typically somewhat smaller. This 
appears to reflect the fragility of the deuteron and its 
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survival probability once formed 18[ . For this reason we 
have used average volumes derived from the A=3 and 4 
clusters to calculate the densities. Given that mass is " 



removed from the system during the evolution, we deter- 1 — Aa'*^ 



mined the relevant masses for each volume by assuming > 

that the initial mass of the source was that determined ^ ■ 



from the source fitting analysis and then determining the I— 
mass remaining at a given Vsmi from the observed energy - ^ - 

spectra. This is also an averaging process and ignores _ ^A.^ 

fluctuations. Densities were determined by dividing re- ^if^ 
maining masses by volumes. 
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Temperatures and Densities 

Inspection of the results for the four different sys- 
tems studied revealed that the temperatures, densities 
for all systems are the same within statistical uncertain- 
ties. Therefore we have combined them to determine the 
values reported in this paper. 

We present, in Figure 1 the experimentally derived 
density and temperature evolution of the IV source. Es- 
timated errors on the temperatures are 10% below p = 
0.01 fm~ increasing to 15% at p = 0.03 fm~3. Estimated 
errors on the densities are 20%. 

In a recently submitted paper we reported equilibrium 
constants for a cluster formation as a function of temper- 
ature and density fl9| . These equilibrium constants were 
then compared with those predicted by several different 
astrophysical equation of state models. Specifically we 
defined the equilibrium constants, Kc, for cluster forma- 
tion in terms of density as 



(3) 



where P(a,z) is the density of clusters of a specific mass 
number A and atomic number Z , N is the neutron num- 
ber in the cluster and pp and p„ are, respectively, the 
densities of free protons and neutrons. In the present 
work, we employ the observed temperature and density 
dependence of these equilibrium constants to extract the 
in medium modifications of the cluster binding energies 
and determine Mott points. 



FIG. 1: Temperatures and densities sampled by the expand- 
ing IV source. 



DERIVATION OF IN MEDIUM BINDING 
ENERGIES 

Our departure point for extraction of the medium mod- 
ified cluster binding energies is the chemical equilibrium 
expression relating the density of a cluster of mass num- 
ber A and atomic number Z to the densities of neutrons 
and protons in the same volume V [l^ . 
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In this expression. At 



(2xmoT)i/^ is the thermal 
wavelength of a nucleon, Sp and s„ are the proton and 
neutron spins, T is the temperature and B{A, Z) is the 
cluster binding energy. The term uj{A, Z) is the internal 
partition function of the cluster, taken here to 1 for the 
Z = 1 and Z = 2 clusters considered. 

Minich et al. [13, [III used a related yield expression 
to analyze intermediate mass fragment yields in multi- 
fragmentation experiments but added both a mixing en- 
tropy term and a surface entropy term, as initially pro- 
posed by Fisher This latter term leads to a power 
law behavior of the mass distribution at the critical point 
and its ramifications have been widely explored [23l425| . 
Neither of these entropy contributions is explicitly in- 
cluded in the Albergo formulation 1^. In the present 
analysis of the experimental data we include a mixing 
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entropy term in the free energy. This term has the form 

l\F = T{Z\n{Z/A) + N\n{N/A)) (5) 

where once again N and A are those of the cluster 
being formed [20| • As mixing is a spontaneous process the 
free energy of mixing is negative and therefore favors the 
cluster formation. We do not include a Fisher term. Our 
reasoning for this is that, without additional corrections, 
the Fisher term as normally formulated and applied to 
larger clusters whose properties are very similar to the 
bulk properties, is not applicable to the yields of the very 
small clusters, A < 4, which we are treating. We base this 
conclusion on the results of molecular dynamics studies of 
the cluster size dependence of the surface energy 26-28j| 
and the binding energies per nucleon of the competing 
Z = 1,2 species being significantly different from the 
bulk. 

Thus, rearranging Eq. (4) , substituting from Eq. 
(3) and taking the logarithm of each side we can write a 
general expression for each cluster, 

\n[KjC{T)] ^B/T~Z \n{Z/A) - N \n{N/A) (6) 

where C{T) includes all terms on the right hand side of 
Eq. Q except the exponential term. Using the experi- 
mentally determined equilibrium constants and tempera- 
tures we then solve this expression to obtain the apparent 
binding energies, B{p, T), of the clusters for the different 
temperatures and densities sampled in the experiments. 
The binding energies extracted for d, t, ^He and a clus- 
ters decrease monotonically with increasing density as 
shown in Fig. [51 



MOTT POINTS 

By definition, a Mott point corresponds to a combi- 
nation of density and temperature at which a cluster 
binding energy, B(p,T), is zero with respect to the sur- 
rounding medium. Since the observed temperatures and 
densities are correlated in our experiment (see Figure 1) 
each point in Figure 2 at which the experimentally de- 
rived binding energy is zero corresponds to a particular 
combination of density and temperature. Thus, with the 
present data, we are able to extract a single Mott point 
for each cluster. In Figure 3 we present the values of 
the Mott temperatures and densities and compare them 
with the loci of the values of medium modified binding 
energies predicted by Typel et al. [H using the thermody- 
namic Green function method. Such a Mott line was also 
calculated in reference [2^ while the contribution of cor- 
relations was considered without discriminating among 
different clusters. This approach makes explicit use of 
an effective nucleon-nucleon interaction to account for 
medium effects on the cluster properties [s^l- We see 
that the agreement between the predictions and the ex- 
perimental results is quite good. 
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FIG. 2: In medium binding energies derived from the exper- 
iments as a function of density. T and p are changing in a 
correlated fashion. (See text.) 
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FIG. 3: Comparison of experimentally derived Mott point 
densities and temperatures with theoretical values. Symbols 
represent the experimental data. Estimated errors on the 
temperatures are 10% and on the densities 20%. Lines show 
polynomial fits to the Mott points presented in reference T|. 
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SUMMARY AND CONCLUSIONS 

We have presented a first experimental determination 
of in medium cluster binding energies and Mott points for 
d, t, '^He and a clusters produced in low density nuclear 
matter. Our results are in good agreement with those 
predicted by a recent model which explicitly treats these 
quantities. Inclusion of the in-medium effects in astro- 
physical equations of state should improve the utility of 
those for modeling astrophysically interesting events. 
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